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Abstract. We present the analysis of a high-resolution, high-quality sample of optical spectra for 76 subdwarf B (sdB) stars 
from the ESO Supernova la Progenitor Survey (SPY, Napiwotzki et al. 2001). Effective temperature, surface gravity, and pho- 
tospheric helium abundance are determined simultaneously by fitting the profiles of hydrogen and helium lines using synthetic 
spectra calculated from LTE and NLTE model atmospheres. We perform a detailed comparison of our measurements with the- 
oretical calculations, both for single star evolution and for binary population synthesis models of close binary evolution. The 
luminosity evolution given by the standard EHB evolutionary tracks from Dorman et al. (1993) shows an overall agreement in 
shape with our observations, although a constant offset in luminosity exists. The various simulation sets for binary formation 
channels of sdB stars calculated by Han et al. (2003) are compared individually to our data for testing our current understanding 
of sdB formation processes and the physical effects involved. The best-matching sets manage to reproduce the observed sdB 
distribution in the temperature-gravity-plane well. However, they do not match the observed cumulative luminosity function, 
indicating that theoretical improvement is necessary. We also investigate composite-spectrum objects showing clear signatures 
of a cool companion with optical and infrared photometry. These stars have cool main sequence companions of spectral types 
F to K. Typical helium abundances of composite and non-composite sdB stars do not differ. 

Key words, binaries: spectroscopic - stars: abundances - stars: atmospheres - stars: fundamental parameters - stars: horizontal 
branch - subdwarfs 



1. Introduction 

Early studies provided evidence that subdwarf B (sdB) stars 
are core helium burning stars with a canonical mass of M x 
0.5 Mq, and a very thin hydrogen envelope (Menv < 0.01 Mq), 
placing them on the very hot end of the horizontal branch (HB), 
the so-called extreme horizontal branch (EHB). Unlike the typ- 
ical post-HB evolution, these objects do not ascend the asymp- 
totic giant branch (AGB) after core helium exhaustion, since 
hydrogen is not burned continuously in a shell due to the very 
low envelope mass. Instead, they evolve more or less directly to 
the white dwarf stage. Proposed formation scenarios are based 
on a late core helium flash in single star evolution (e.g. D'Cruz 
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et al. 1996) or on mass transfer in close binary systems (e.g. 
Mengel et al. 1976); furthermore, the merger of two helium 
white dwarfs (He-WDs) could be the origin of sdB stars (Iben 
1990). 

Today, more than 200 sdB stars have been analyzed for at- 
mospheric parameters (e.g. Saffer et al. 1994; Maxted et al. 
2001; Edelmann et al. 2003), and this interpretation still seems 
to be valid in general: effective temperature and surface grav- 
ity of the observed sdBs mostly lie within the theoretically 
determined start and end points for core helium burning on 
the EHB, zero-age EHB (ZAEHB) and terminal-age EHB 
(TAEHB). However, the extensive theoretical study of binary 
formation mechanisms of sdB stars by Han et al. (2003, here- 
after HPMM) shows that a huge number of sdB stars may be 
missing in observational surveys due to selection effects. These 
are primarily caused by main sequence companions that out- 
shine the sdBs (main sequence spectral type A and earlier) or 
appear as composite spectrum objects (main sequence type F to 
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K). The latter have mostly been set aside because their spectral 
analysis is rendered difficult by the companion spectrum. 

More than ten years after the first quantitative estimate of 
the contribution of different binary channels to the population 
of sdB stars (Tutukov & Yungelson 1990), it is clear from the 
observations of Maxted et al. (2001) that close binary evolu- 
tion is indeed of great importance to sdB formation processes, 
since they find that two thirds of their observed sdB stars are 
members of close binaries. The variety of physical parameters 
involved in those processes led HPMM to produce twelve sim- 
ulation sets for sdB stars that formed in several evolutionary 
channels, i.e. by stable Roche lobe overflow (RLOF), com- 
mon envelope (CE) ejection, or merging of two He-WDs. In 
comparing them with an observational sample providing good 
quality in atmospheric parameters, it should be possible to con- 
strain the simulations' parameter range and to draw conclu- 
sions about the relative importance of each formation channel. 

The ESO Supernova la Progenitor Survey (SPY, 
Napiwotzki et al. 2001) provides such a sample, since 
the ambitious project obtained high-resolution optical spectra 
of over 1000 white dwarf candidates, containing some 140 
previously misclassified hot subdwarfs of various types. Earlier 
analyses of sdBs (e.g. Saffer et al. 1994; Maxted et al. 2001; 
Edelmann et al. 2003) were limited by the lower resolution of 
their spectra, sometimes drawn from inhomogeneous data sets. 
In many cases the wavelength coverage was incomplete as 
well, e.g. not all Balmer hues were included. Here we present 
the quantitative analysis of an unprecedented homogeneous set 
of high-resolution, high-quality spectra with large wavelength 
coverage, all of which makes it an excellent means for testing 
and increasing the knowledge of sdB formation and evolution. 

Our report on the analysis of 76 sdB stars from the SPY 
sample is outlined as follows. In Sect. 2, we describe the ob- 
servations and data reduction of our objects. Section 3 contains 
information on the model atmospheres and the fit procedure 
used to determine atmospheric parameters, as well as on the 
examination of the Ha line and of spectral Unes from cool com- 
panions. In Sect. 4, we analyze the sdB stars showing no spec- 
tral signatures of a cool companion and compare them to ear- 
Uer studies. The composite-spectrum objects are investigated 
in Sect. 5, where we attempt to draw conclusions about the na- 
ture and spectral types of the companions. Finally, in Sect. 6, 
a detailed comparison of our observations with calculations for 
sdB formation and evolution is presented. 

2. Observations and data reduction 

Observations were obtained at the ESO Very Large Telescope 
with UT2 (Kueyen) equipped with the UV- Visual Echelle 
Spectrograph (UVES, Dekker et al. 2000). A slit width of 2'.'1 
was used, resulting in a spectral resolution of 18 500 (0.36 A 
at Hot) or better. Wavelength coverage of 3300 - 6650 A is 
achieved, with gaps at 4500 - 4600 A and 5600 - 5700 A. For 
most of the stars, two exposures in different nights were taken, 
since SPY was originally intended to search for RV-variable 
objects. The spectra were then reduced with a procedure de- 
veloped by Karl (in prep.) using the ESO MIDAS software 
package, partly based on the UVES onhne reduction facility. 



After accounting for cosmic ray hits and bad CCD pixels, bias 
and interorder background were subtracted, followed by order 
extraction for sky background, flatfield, and object. For each 
extracted one dimensional order, sky background subtraction 
and flatfielding were performed, and finally the orders were 
merged, resulting in three partial spectra separated by the gaps 
mentioned above. Those spectral parts are hereafter referred 
to as blue, lower red, and upper red part, respectively. Finally, 
each of them was divided by a smoothed spectrum of a DC 
white dwarf, which by definition shows no spectral features at 
all and therefore provides an excellent means of correcting for 
the instrumental response. 

The spectra were then convolved with a Gaussian of 1.0 A 
FWHM and rebinned to 0.4 A stepsize. A signal-to-noise ra- 
tio (S /N) was defined and calculated for the blue parts of the 
spectra, since they contain most of the Unes used for fitting, 
as the ratio of mean flux to standard deviation in a continuum 
area. The median value is 178, all except one have values of 
> 10, nine spectra even have S/N> 1000. 

Photometric data are available for all of our objects. B mag- 
nitudes of most of the stars are photographic measurements 
provided by the Hamburg/ESO Survey (HES, Wisotzki et al. 
1996), the Hamburg Quasar Survey (HQS, Hagen et al. 1995), 
and McCook & Sion (1999). For nine HES stars, and V val- 
ues were taken from the CCD photometry of Altmann et al. 
(2004). In some cases, only V or Stromgren-y magnitudes were 
available. J magnitudes have been drawn from the Two Micron 
All Sky Survey (2MASS) archive. 

3. Spectral analysis and classification 

3.1. Model atmospheres and line profile fitting 

Effective temperatures (Tetf), surface gravities (log (g)), and he- 
hum abundances (y = Nhc/Nh) were determined by fitting 
simultaneously each hydrogen and helium line to synthetic 
model spectra, using a procedure developed by R. Napiwotzki 
(Napiwotzki et al. 1999) based on Saffer et al. (1994). Three 
different sets of models were used, where the last two make 
use of PR02 (Werner & Dreizler 1999): 

1. A grid of metal line-blanketed LTE model atmospheres 
with solar abundance (Heber et al. 2000) for stars with 
Teff < 32 000 K. 

2. A grid of partially line -blanketed NLTE model atmospheres 
(Napiwotzki 1999) for stars with T^b > 32 000 K. 

3 . A grid of partially line-blanketed NLTE model atmospheres 
for helium rich objects, used for two stars with log (y) > - 1 
(updated version of Dreizler et al. 1990). 

A typical LTE and NLTE fit is shown in Fig. 1. For each 
hne, the continuum level is determined and normalized to 1, in 
order to compare it to the synthetic spectrum. Ha was never 
included in the parameter determination itself, but was kept in 
the final plot for examining possible deviations from the model 
Ha line. The core of Hor is sensitive to NLTE effects because it 
is formed in the very outer parts of the atmosphere. Therefore, 
LTE models often fail to reproduce the observations. This is ev- 
ident in the example given for a LTE analysis in Fig. la. NLTE 
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Fig. 1. a) Sample fit of a relatively cool sdB (GD 687) using synthetic 
spectra calculated from LTE model atmospheres, b) Sample fit of a 
relatively hot sdB (GD 619) using synthetic spectra calculated from 
NLTE model atmospheres. In both panels, the observed spectra are 
shown as histograms. When a spectral line was used for parameter 
determination, model spectra are plotted as solid grey lines, otherwise 
as dashed grey lines. The continua of observed and model spectrum 
were scaled to the same level on both sides of the line. The resulting 
atmospheric parameter values are displayed in the figure. 



analyses, however, usually are able to reproduce the Ha pro- 
files of sdB stars very well (Heber et al. 2003), as exemplified 
in Fig. lb. We will discuss the Ha lines of the programme stars 
in more detail in Sect. 3.4. 



Fig. 2. Colour-magnitude diagram of B versus B - J for all 76 ob- 
jects. Composite spectrum objects are plotted as grey circles if they 
are classified from the presence of the Mg i triplet in their spectra, or 
as grey diamonds if they are classified from the contribution of the 
companion to the Ha line profile (see text). Asterisks represent the 
four stars with peculiar Ha profiles. The dotted vertical line marks the 
value B - J - (FO. For the stars marked as open symbols, we adopted 
B X V - 0P25 or fi K y - 0?25, because no B measurements are 
available. 



3.2. Spectral classification 

The hot subdwarfs had been selected from the SPY data by us- 
ing Balmer line widths as an indicator for surface gravity and 
by performing a line profile fit where visual inspection did not 
yield unambigous results. After having excluded objects that 
clearly have helium-dominated atmospheres, the presence and 
strength of the He n lines was used for distinguishing between 
sdB, sdOB, and sdO stars. In this paper we present analyses of 
the sdB and sdOB stars from SPY, subsequently referred to as 
sdB stars. This sample comprises 76 sdBs with effective tem- 
peratures between 20 000 K < Tetf < 38 000 K, gravities in the 
range 4.8 < log (g) < 6.0, and helium abundances between 
-4.0 < \og(y) < -0.8 (see Table 1). 51 program stars are ob- 
jects from HES, 9 from HQS, 14 from McCook & Sion (1999), 
and 2 are from the hst of white dwarf central stars of planetary 
nebulae by Napiwotzki (1999), see Table 1. 

3.3. Spectroscopic and photometric signatures of cool 
companions 

Our sample of sdB stars includes many objects where signa- 
tures of a cool companion can be seen in the spectrum. In 
some of them, the companion shows a large contribution to 
the total flux even in the blue part of the spectrum, causing 
a broad Ca K line (3933 A) as opposed to the frequently occur- 
ring narrow interstellar Ca K Une. In addition, G-band absorp- 
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Table 1. Results of our spectral analysis. T^g, log (g) and log (y) are mean values from two or more exposures, if possible. Luminosity in units 
of solar luminosity is calculated by assuming = 0.5 M©. 



Object 




log (e) 


log (v) 


log(L) 


log(L) 


My 


d 


RA(J2000) 


DEC (J2000) 


B 


Note 




K 


cm s~^ 








mag 


kpc 


hh:mm:ss.ss 


±dd:mm:ss.s 


mag 




EGB5 


34060 


5.85 


-2.77 


-2.84 


1.37 


4.50 


0.7 


08:11:12.77 


+ 10:57:16.8 




7 


HE 0007-2212 


28964 


5.68 


< -4.00 


-2.95 


1.25 


4.42 


1.4 


00:09:45.91 


-21:56:14.4 


15.0 


10 


HE 0016+0044 


28264 


5.38 


-2.66 


-2.70 


1.51 


3.72 


1.7 


00:18:43.59 


+01:01:22.5 


14.8 




HE 0019-5545 


35662 


5.86 


-1.44 


-2.76 


1.44 


4.45 


2.2 


00:21:27.68 


-55:29:12.3 


15.8 




HE 0101-2707 


35568 


5.97 


-0.92 


-2.88 


1.32 


4.73 


1.2 


01:03:43.73 


-26:51:53.8 


15.1 


10 


HE 0123-3330 


36602 


5.87 


-1.49 


-2.73 


1.47 


4.45 


1.6 


01:25:22.46 


-33:15:11.4 


15.2 


10 


HE 0135-6150 


27020 


5.59 


-2.47 


-2.99 


1.22 


4.34 


2.8 


01:37:30.08 


-61:34:57.8 


16.3 




HE 0136-2758 


28202 


5.47 


< -4.00 


-2.79 


1.41 


3.95 


2.8 


01:39:14.46 


-27:43:21.7 


15.94 


5,10 


HE 0151-3919 


20841 


4.83 


-2.07 


-2.68 


1.53 


2.94 


1.9 


01:53:11.18 


-39:04:18.2 


14.09 


5,9 


HE 0207+0030 


31414 


5.83 


-1.83 


-2.95 


1.25 


4.63 


1.4 


02:10:14.97 


+00:45:02.1 


15.1 


8 


HE 0230-4323 


31552 


5.60 


-2.58 


-2.72 


1.48 


4.04 


0.9 


02:32:54.66 


-43:10:27.9 


13.56 


5 


HE 0306-0309 


26710 


5.79 


< -4.00 


-3.20 


1.01 


4.86 


1.9 


03:08:40.80 


-02:58:02.8 


16.1 




HE 0321-0918 


25114 


5.67 


-3.02 


-3.18 


1.02 


4.69 


1.3 


03:23:45.85 


-09:08:15.6 


15.2 




HE 0415-2417 


32768 


5.12 


-2.44 


-2.18 


2.03 


2.76 


5.6 


04:17:31.49 


-24:09:50.7 


16.2 




HE 0513-2354 


26758 


5.50 


-2.33 


-2.91 


1.29 


4.14 


2.4 


05:15:15.30 


-23:51:09.5 


15.8 




HE 0532-4503 


25710 


5.33 


-3.07 


-2.81 


1.40 


3.79 


2.8 


05:33:40.51 


-45:01:35.3 


15.84 


5 


HE 0539-4246 


23279 


5.51 


-3.91 


-3.16 


1.05 


4.43 


1.2 


05:41:06.71 


-42:45:31.9 


14.60 


5 


HE 0929-0424 


29602 


5.69 


-2.01 


-2.92 


1.28 


4.40 


1.9 


09:32:02.15 


-04:37:37.8 


15.4 




HE 1021-0255 


35494 


5.81 


-1.46 


-2.73 


1.48 


4.33 


1.7 


10:24:18.06 


-03:10:38.8 


15.3 


1,8 


HE 1033-2353 


36204 


5.76 


-1.46 


-2.64 


1.56 


4.19 


2.6 


10:36:07.23 


-24:08:35.4 


16.0 


9 


HE 1038-2326 


30573" 


5.21'' 


-1.73" 










10:40:36.97 


-23:42:39.4 


15.9 


1,2,9 


HE 1047-0436 


30280 


5.71 


-2.35 


-2.91 


1.30 


4.41 


1.3 


10:50:26.93 


-04:52:35.8 


14.7 


8 


HE 1050-0630 


34501 


5.79 


-1.40 


-2.76 


1.45 


4.33 


0.9 


10:53:26.52 


-06:46:15.6 


14.0 


8 


HE 1 140-0500 


34522" 


4.97" 


-2.62" 










11:42:57.85 


-05:17:14.1 


14.8 


2,8 


HE 1200-0931 


33419 


5.78 


-1.88 










12:03:21.77 


-09:48:06.6 


16.2 


3 


HE 1221-2618 


32606" 


5.51" 


-1.84" 










12:24:32.72 


-26:35:17.1 


14.7 


2,9 


HE 1254-1540 


29700 


5.63 


< -4.00 










12:57:19.36 


-15:56:22.8 


15.2 


3,9 


HE 1309-1102 


27109" 


5.36" 


-2.46" 










13:12:02.39 


-11:18:16.3 


16.1 


2 


HE 1352-1827 


35674" 


5.53" 


-1.77" 










13:55:26.67 


-18:42:09.4 


16.0 


2,9 


HE 1407+0033 


37309 


5.54 


-2.99 


-2.38 


1.83 


3.61 


2.4 


14:10:20.73 


+00:18:54.6 


15.5 


4,8 


HE 1415-0309 


29520 


5.56 


-2.87 


-2.80 


1.41 


4.09 


3.0 


14:18:20.93 


-03:22:54.1 


16.3 




HE 1419-1205 


34171 


5.71 


-1.65 










14:22:02.17 


-12:19:30.9 


16.2 


3 


HE 1421-1206 


29570 


5.53 


< -4.00 


-2.76 


1.45 


4.01 


1.9 


14:24:08.81 


-12:20:21.5 


15.1 


8 


HE 1422-1851 


33896" 


5.19" 


-3.07" 










14:24:48.65 


-19:05:01.3 


16.3 


1,2,9 


HE 1441-0558 


36396" 


5.79" 


-1.63" 










14:44:12.11 


-06:10:44.7 


14.4 


2 


HE 1448-0510 


34760 


5.53 


-3.41 


-2.48 


1.72 


3.66 


1.5 


14:51:13.13 


-05:23:16.9 


14.4 


4,8 


HE 1450-0957 


34563 


5.79 


-1.29 


-2.75 


1.45 


4.32 


1.7 


14:53:24.19 


-10:09:21.9 


15.1 


9 


HE 1459-0234 
















15:02:12.29 


-02:46:00.9 


14.9 


2 


HE 1519-0708 


34498 


5.73 


-1.52 


-2.70 


1.50 


4.18 


2.3 


15:21:53.20 


-07:19:23.6 


15.6 


8 


HE 2135-3749 


29924 


5.87 


-2.45 


-3.08 


1.13 


4.83 


0.7 


21:38:44.18 


-37:36:15.1 


13.70 


5 


HE 2150-0238 


29846 


5.90 


-2.36 


-3.12 


1.09 


4.91 


1.7 


21:52:35.81 


-02:24:31.6 


15.8 




HE 2151-1001 


34984 


5.70 


-1.60 


-2.64 


1.57 


4.07 


2.2 


21:54:31.49 


-09:47:30.5 


15.6 


4 


HE 2156-3927 


27995" 


5.50" 


-2.35" 










21:59:35.53 


-39:13:15.3 


14.26 


2,5,11 


HE 2201-0001 


27062 


5.51 


-3.29 


-2.90 


1.31 


4.14 


2.7 


22:04:18.27 


+00:12:36.7 


16.0 




HE 2208+0126 


24277 


5.67 


-2.98 


-3.25 


0.96 


4.75 


1.4 


22:10:45.47 


+01:41:35.4 


15.2 


1,8 


HE 2222-3738 


30248 


5.69 


-3.65 


-2.88 


1.32 


4.36 


1.3 


22:24:56.50 


-37:23:30.7 


14.66 


5 


HE 2237+0150 


25606 


5.38 


-1.92 


-2.86 


1.34 


3.92 


2.7 


22:40:14.38 


+02:06:31.3 


15.8 


8 


HE 2238-1455 


30393 


5.47 


-2.37 


-2.66 


1.54 


3.80 


3.0 


22:41:38.27 


-14:39:39.5 


16.0 




HE 2307-0340 


23260 


5.51 


-3.65 


-3.16 


1.04 


4.44 


2.1 


23:10:24.09 


-03:24:02.3 


15.8 




HE 2322-0617 


28106" 


5.50" 


-1.93" 










23:25:31.91 


-06:01:12.1 


15.7 


2 


HE 2322-4559 


25512" 


5.30" 


-2.47" 










23:25:09.05 


-45:43:06.5 


15.5 


2 


HE 2349-3135 


28520 


5.44 


-3.84 


-2.74 


1.47 


3.86 


2.6 


23:51:43.63 


-31:18:52.9 


15.63 


5,10 


HS 1530+0542 
















15:33:10.74 


+05:32:26.8 


14.3 


2,8 


HS 1536+0944 


35114" 


5.83" 


-0.82" 










15:38:42.88 


+09:34:42.8 


15.6 


1,2,8 
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Object 


7'eff 


log (g) 


logO) 


log(L) 


log(L) 


Mv 


d 


RA (J2000) 


DEC (J2000) 


B 


Note 




K 








Lo 


mag 


kpc 


hh:mm:ss.ss 


±dd:mm:ss.s 


mag 




HS 1710+1614 


34826 


5.72 


-1.64 


-2.68 


1.53 


4.13 




i / . 1 j.Uj. lo 


-1-1 A-in-zt9 s 

T lO. ivJ.HZ.O 


1 ^ 7 
1 J. / 


A 


HS 2033+0821 


32706 


5.87 


-1.56 


-2.93 


1.28 


4.64 


i.U 


zu. jj.zy. j^- 


+UO.J1.J1. / 


^AA 




HS 2043+0615 


26157 


5.28 


-2.38 


-2.73 


1.48 


3.63 




ZU.T-D.ZU.OD 


j-nfi.Ofi.O/i A 

T^UU.ZO.ZH-.H- 


1 n 
10. U 




HS 2125+1105 


32542 


5.76 


-1.86 










Zl .z / . JZ. 1 / 


-.11.10.17 1 

Tll.lO.l / .I 


ID.T- 


-J 
J 


HS 2216+1833 


3436P 


5.51'' 


-1.70" 










ZZ. i o. 3U.UV^ 


i O.HO.UV.H 


LJ.O 


9 
z 


HS 2357+2201 


27629 


5.55 


-2.54 


-2.90 


1.30 


4.20 
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" The given value is an upper limit due to the presence of a cool companion. 

" The given value is a lower limit due to the presence of a cool companion. 

° Johnson V magnitude 

" Stromgrcn }' magnitude 

1 = stars for which there was only one useful exposure. 

2 = a cool companion shows Mg i in the spectrum, and possibly additional features. 

3 = the presence of a companion was deduced solely from a flux contribution at He Atmospheric parameters have been determined from 
Balmer and helium hnes in the blue part of the spectrum. Unlike for the stars with earlier type companions, the contribution of the 
companion in this spectral region is irrelevant, and the parameters can be regarded as reUable (for details see Sect. 5.1). 

4 = single-lined objects showing pecuhar Hor profiles. 

5 = independently classified as sdB by Altmann et al. (2004). 

6 = central star of the planetary nebula PN G 125.9-47.0 . 

7 = central star of the planetary nebula PN G 211.9+22.6 . 

8 = also in the Palomar-Green survey (PG, Green et al. 1986). 

9 = also in the Edinburgh-Cape survey (EC, Stobie et al. 1987). 

10 = also in the Montreal-Cambridge-Tololo survey (MCT, Demers et al. 1987). 

11 = In Table 1 of Altmann et al. (2004), V magnitudes of this star and HE 2156-1732 have to be interchanged (M. Altmann, priv. comm.). 



tion was taken as a clear companion signature. For most ob- 
jects though, the contribution of the companion is best visible 
in the lower red part of the spectrum, where the Mg i triplet at 
5167 A, 5173 A, and 5184 A is the most striking feature. Since 
the latter is not known to occur in sdB atmospheres nor in the 
interstellar medium, we used it as a tracer to find even very faint 
companions in sdB spectra that appeared to be non-composite 
at first glance. These examinations yielded altogether 19 sdB 
stars with a cool companion contributing to the spectrum (see 
Table 1). They are treated separately (cf. Sect. 5) from the non- 
composite objects, because the results from fine profile fitting 
cannot be considered reUable. 



Since our spectra have not been flux calibrated, it was 
not possible to verify these findings by making use of the 
flux distribution in the upper red part, which should otherwise 
clearly show the contribution of the companion. Instead, we re- 
trieved infrared photometric data for our stars from the 2MASS 
database to examine the quality of the Mg i triplet as tracer for 
cool companions. 67 of our objects show a detection at least in 
J and are unambiguously identified, whereas 9 show no detec- 
tion at all and cannot be mistaken with nearby sources. In the 
latter cases, we derived J < 16!^5 as a conservative estimate of 
the detection Umit from the fluxes of nearby sources in the field 
of view. 
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Fig. 3. Comparison of synthetic Hor line profiles calculated from NLTE model atmospheres to observed ones (see text), a) GD 619: a typical 
sdB, b) HE 1441-0558: a composite sdB star, c) HE 1441-0558: after correction for continuum light from the cool companion. 




-20 

relative wavelength 



Fig. 4. Same as Fig. 3, but here we show the second kind of Ha profile 
peculiarity we observed, in the spectrum of HE 1448-0510 (shown 
here) and four other stars (for details see text). 



The 2MASS observations are complete^ down to 7 < 15?8, 
H < 15^1 and < 14?3. Since more than half of our sample is 
fainter than J - 15"5, their H and K fluxes are of little use and 
were not considered. We based our analysis onB - J colours, 
because V-band measurements are unavailable for most stars. 
In Fig. 2 we plot Bvs.B-J for all 76 sdB stars. The apparently 
single sdB stars separate from objects with a cool companion 
at about B - 7 « 0, which nicely agrees with the detection of 
Mgi in their spectra (grey circles in Fig. 2). 

For companion spectral types later than KO, the differ- 
ence in colours between an sdB+main sequence system and 
an sdB alone decreases and finally almost vanishes, because of 
the decreasing contribution of the companion to the total flux 
(cf. Sect. 5.1). Hence, there is a limit in spectral type above 
which an sdB+main sequence system may lie at B - / < 0, and 
therefore cannot be identified as composite. Three stars with 
B - J < show Mgi as weU, confirming the quality of the 
latter as an even better tracer for cool companions than colour 
values are. 



' http://www.ipac.caltech.edu/2mass/overview/about2mass.html 



3.4. Peculiar Ha line profiles 

A close examination of the Ho- fine of the apparently non- 
composite objects yielded ten stars with pecuUar Ha profiles. 
They can be divided into two groups: five of them show a Une 
core which is shallower than that of the corresponding model 
Une, but has a normal fine core shape. This peculiarity is seen 
in seven composite objects as well, suggesting that this is an- 
other hint of a cool companion. It is demonstrated in Fig. 3 
where we compare synthetic Ha Une profiles calculated from 
NLTE model atmospheres to observed ones for a typical sdB 
(Fig. 3a) and for a representative composite one showing this 
peculiarity (Fig. 3b). Atmospheric parameters have been de- 
rived for each star from other B aimer and helium lines in the 
blue part of its spectrum. It must be stressed that we did not fit 
the Ha Une profile by a synthetic spectrum. As can be clearly 
seen, the Ha Une profile in a typical sdB is quite well matched, 
whereas in the pecuUar star it is too shallow. Since the cool 
companion adds some continuum flux, it dilutes the Ha pro- 
file of the sdB. Hence, after normalization, the Ha line depth 
is lower than it would be without the companion contribution, 
and therefore lower than that of the synthetic profile. This is 
exactly what we see in those stars. We subtracted an appropri- 
ate constant from the spectrum around Ha. Now the corrected 
Ha profile is perfectly matched by the synthetic profile (see 
Fig. 3c). Since this is shown for a star that we already classi- 
fied as composite, by subtracting a constant we automatically 
demonstrated that the contribution of the companion Ha Une is 
negligible. Applying this treatment to the five non-composite 
stars leads to the same result, which we take as evidence for 
the presence of a cool companion (see Table 1), although other 
spectral features (i.e. Mg i) are not detected in these stars. 

Since the companion contribution decreases with decreas- 
ing wavelength, it is understandable why in the blue and lower 
red part of the spectrum no cool star features can be seen. The 
companions are therefore expected to be of very late type, con- 
sistent with the colour-magnitude-diagram (B, B - J) (Fig. 2): 
one star lies among the composite objects, the others lie among 
the single sdBs, but are very close to B - 7 = 0, which we re- 
gard as the photometric detection limit for composite objects. 
As discussed in the previous subsection, this is exactly what 
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Table 2. Errors derived from the Gaussians in Fig. 5. For stars with 
one exposure (last column), the errors equal the Gaussian widths. For 
objects with two exposures or more, those values are divided by V2 
(middle column). 



Quantity 


Errors 


Errors " 


T^ff 


374 K 


529 K 


log (g) 


0.049 dex 


0.069 dex 


logCv) 


0.044 dex 


0.062 dex 


logCL)" 


0.038 dex 


0.053 dex 



" For stars with only one useful exposure. 
* Luminosity L. 



one expects for very late type companions, thus supporting our 
conclusion. We summarize that 24 out of 76 sdB stars show 
signatures of a cool companion, hereafter termed double-lined 
objects, even if the stars just discussed do not show any ob- 
vious companion lines. As mentioned above, they are treated 
separately from the 52 non-composite (single-lined) sdBs. 

A second kind of peculiarity of the Ha profile is found in 
five other objects of our sample: the core is flat, and the in- 
ner wing is broader than predicted by the model, as seen in 
Fig. 4. One of these stars, EGB 5, is the central star of a plane- 
tary nebula. The Ha emission of its nebula may not have been 
correctly subtracted in the semi-automatic data reduction pro- 
cess. Close inspection of sdB spectra that we regard as normal 
(Fig. 3a) reveals that slight mismatches of the Ha line core oc- 
cur even in those stars. However, these effects are so small that 
they could be caused by errors in the data reduction and/or the 
atmospheric models and are therefore considered as insignifi- 
cant. Nevertheless, we are left with four stars with peculiar Ha 
profiles. 

The colour values of this group of pecufiar stars (Fig. 2, 
asterisks) clearly show that they cannot be composite objects. 
There is no significant line profile variation observed from one 
to the next exposure of the star. Time intervals range from 3 d 
to 1 yr. Although the observed line broadening could possibly 
be caused by rotation, convolution of the synthetic spectra with 
rotational profiles does not lead to acceptable fits for any hy- 
drogen or helium line. Even for Ha itself, the match is still not 
satisfying, implying that other physical effects are present. 

Since at this point there is no interpretation for this Ha pro- 
file peculiarity, we proceed with our analysis, and will distin- 
guish those stars from the rest of single-lined objects wherever 
appUcable. 

4. Single-lined objects 

4.1. Error determination 

The statistical l-cr-errors from the fit procedure are typically 
lower than lOOK, 0.02 dex, and 0.04 dex for T^g, logig), and 
log (y), respectively. However, despite the high resolution and 
the low noise level of our data, we do not believe these to be 
the "true errors". To obtain a good estimate for the latter, we 
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Fig. 5. Distributions of differences in atmospheric parameters between 
two exposures of non-composite stars. Solid lines show the Gaussians 
fitted to the data. From top to bottom the distributions for Tgg, log (g), 
log (y) and luminosity L with respect to the Eddington luminosity are 
shown. 

take advantage of the SPY observing strategy: as mentioned 
above, at least two exposures were taken for most of the stars. 
Hence, for determining reliable errors we used the distribution 
of differences in the fit results of individual exposures of each 
object (Fig. 5). Fitting this distribution by a Gaussian, we adopt 
the width of the latter as the error of the fit results of individ- 
ual exposures. This is a quite conservative treatment intended 
to avoid underestimating the errors: since it is usually assumed 
that the individual results show a Gaussian error distribution as 
well, the width of the latter would be lower by a factor of y/l 
than the width we derived from our distribution of differences 
of the results of two exposures. However, we do not think that 
this standard assumption would be well justified here, since in 
many cases the differences are much larger than the typical sta- 
tistical errors mentioned above. 

Furthermore, when we take the mean value of two expo- 
sures, the errors are automatically reduced by a factor of V2, 
increasing the accuracy of the derived parameters. For stars 
with two or more spectra, those reduced Gaussian widths are 
ATf.f[ - 360 K and A log (g) - 0.05 dex, which we regard as 
our best error estimates. More details are given in Table 2. 

An uncertainty is caused by the unknown metallicity of the 
program stars. We used solar metallicity, fully line-blanketed 
LTE models (see Sect. 3.1). A comparison shows that adopt- 
ing a much lower metallicity ([m/H] = -2) in the LTE models 
has a negligible effect at T^g - 32 000 K (see also Heber et al. 
2000) but increases T^g by up to 800 K at lower temperatures. 
Surface gravity is only marginally affected (0.02 dex) while the 
helium abundance remains almost unchanged. Spectroscopic 
analyses of high resolution optical spectra of about two dozen 
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Fig. 6. Distribution of our sdB sample in the Teif-log (g)-plane. Results 
derived from only one exposure are shown as open circles, whereas 
filled circles represent mean values from two or more spectra. Errors 
are quoted separately for the former and the latter. Asterisks represent 
the four stars with peculiar Ho- profiles. The ZAEHB (short-dashed 
line), TAEHB (long-dashed line) and He-ZAMS (dashed-dotted line) 
are shown (see text for details). 



bright sdB stars found near-solar iron abundances (Heber et al. 
2000; Heber & Edelmann 2004), corroborating our choice of 
metallicity. Hence the uncertainties introduced by the unknown 
metallicity are smaller than the errors of the spectroscopic anal- 
yses quoted above for the majority of our objects. 

Another source of systematic error lies with the usage of 
both LTE and NLTE model atmospheres. Heber & Edelmann 
(2004) have shown that a systematic offset in gravity exists: 
analysis with our grid of NLTE models leads to a gravity lower 
by 0.06 dex than derived from our LTE models, while the re- 
sulting effective temperatures and helium abundances are the 
same. 

We find that neither of the effects above affects any of our 
conclusions (see Sect. 6). 

4.2. Effective temperature and surface gravity 

The fit results are listed in Table 1 . They were calculated as 
the mean value of the results from individual exposures of each 
star, using only spectra with S /N > 10. For nine stars, only one 
(useful) exposure was available (see Table 1). 

Figure 6 shows the position of our sdB stars in the Teff- 
log (^)-plane, along with the ZAEHB and TAEHB for solar 
metallicity (Dorman et al. 1993) and the theoretical zero-age 
main sequence for pure helium stars (He-ZAMS, Paczynski 
1971). Up to effective temperatures of Teff ~ 32 000 K, most of 
our sdB stars lie in the so-called EHB strip, defined by ZAEHB 
and TAEHB. For higher Teg, almost all of the stars lie near the 




log(y) 

Fig. 7. Shaded histogram: Distribution of photospheric helium abun- 
dance of our sdB stars. Open histogram: Distribution from Edelmann 
et al. (2003). The value -4 is the lowest helium abundance in the 
model atmospheres of both studies. Since no extrapolation was done, 
it has to be regarded as an upper Umit for the true abundance. 

TAEHB or above it, including the four objects with peculiar 
Ha profiles. Six stars fall slightly below the ZAEHB. Possible 
explanations could be selection effects and/or evolutionary ef- 
fects, which will be examined in Sects. 4.5 and 6, respectively. 

From T^g and log (g), the luminosity could be calculated 
when a mass for the sdB is assumed. However, we prefer to use 
luminosity in units of the Eddington luminosity, thereby elimi- 
nating the mass. It is calculated from 

log(j^) = 4 X log(Teff/K) - log(^/cms-2) - 15.118 , 
where pure electron scattering in a fuUy ionized hydrogen at- 
mosphere is assumed. 

In order to estimate distances we derive the absolute vi- 
sual magnitude My from the surface size and flux. To calculate 
the former, we need the surface gravity, and an assumption for 
the mass, taken to be MgdB = 0.5 Mq. The surface flux is cal- 
culated from Kurucz model atmospheres (Napiwotzki 2001). 
Subsequently, the distance is derived from My and the appar- 
ent visual brightness my. For most stars, though, only /wb is 
available. Therefore, we adopted B - V = -0^25 for these ob- 
jects, which is the typical value for non-composite sdB stars 
(Stark & Wade 2003, see also Fig. 2). We assumed Stromgren- 
y to be identical to V. All these additional quantities are given 
in Table 1. 

4.3. Helium abundance 

Our distribution of helium abundances (Fig. 7) shows a wide 
spread from the lower boundary of our model grids, log (y) - 
-4, with sometimes no helium visible at all, to slightly super- 
solar helium abundance log (y) = -0.92. There is overall agree- 
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Fig. 8. Helium abundance vs. elfective temperature. Asterisks repre- 
sent the four stars with peculiar Hcc profiles. Arrows denote upper 
liiruts for helium abundance. The dotted horizontal line marks a value 
of log (y) = -2, the vertical line Teff = 32 000 A". The dashed Unes are 
from Edelmann et al. (2003), see text for details. 
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Fig. 9. Helium abundance vs. luminosity with respect to the Eddington 
luminosity. Asterisks represent the four stars with peculiar Ha pro- 
files. Arrows denote upper Umits for helium abundance. The inset 
shows the luminosity distribution as histogram for the four stars with 
peculiar Hof profiles (shaded) and the other objects (open). 



ment with the results from the sample of Edelmann et al. (2003) 
(see Fig. 7), but we find a higher fraction of stars with very low 
helium abundances. Possible correlations with other quantities 
are presented in Figs. 8 and 9, where helium abundance is com- 
pared to effective temperature and luminosity, respectively. The 
comparison with Teff shows that for log (y) > -2, almost all of 
the sdB stars are hotter than Teff > 32 000 K (see the dotted 
lines in Fig. 8). At these high temperatures, only two objects 
have a very low helium abundance \o%{y) < -3.0, whereas 
there are twelve such stars at lower temperatures. Hence, at 
TeB ~ 32 000 K, our objects could be divided into two groups 
that are roughly correlated with spectral type sdB and sdOB, 
respectively, the latter having higher temperatures and an over- 
all higher helium abundance. 

In addition, a slight trend with luminosity can also be ob- 
served: At higher Teff, the sdB stars are somewhat more lu- 
minous (see Table 1). Consequently, the average luminosity is 
higher for higher helium abundances, as seen in Fig. 9. 

Interestingly, the two stars which have the lowest helium 
abundance of the high-temperature objects show peculiar Ha 
profiles, and are more luminous than the other two stars of this 
kind. One might thus speculate about underlying physical ex- 
planations. However, the Ha profiles of these two subgroups do 
not show any obvious differences, and from a statistical point 
of view there would be no significance in such a separation. 

Edelmann et al. (2003) found two distinct sequences of 
sdB stars, separated mainly by an offset in helium abundance. 
Their corresponding linear fits are also plotted in Fig. 8. The 
sequence for higher helium abundance (short dashed line) 
matches the majority of our objects quite well, especially 



at higher temperatures. The sequence for lower log (y) (long 
dashed line), with which Edelmann et al. (2003) matched about 
l/6th of their data, can be neither confirmed nor rejected, 
mainly because we can only derive upper limits for five stars. 

As discussed in previous work (e.g. Fontaine & Chayer 
1997), the diffusive equilibrium abundances of gravitational 
settling versus radiative levitation lie two orders of magnitude 
below the observed average helium abundances; hence, this 
diffusion model is too simplistic. Fontaine & Chayer (1997) 
included mass loss of IQ-^^MQ/yv as an additional force in 
their diffusion models, and concluded that it is possible to reach 
equilibrium abundances as high as log (y) - -3 after t - 10^ yr 
on the EHB. Unglaub & Bues (2001) showed that mass loss 
rates of lO^'^MQ/yr can produce even higher helium abun- 
dances log (y) = -2, again after t - 10^ yr. The large number 
of stars with log (y) > -2 in our sample indicates that either the 
winds must be even stronger for reproducing the whole range 
of abundances, or alternative explanations have to be found. 

Since the stellar winds of sdB stars are probably radiation 
driven, their mass loss rates are expected to be correlated with 
their luminosities (Pauldrach et al. 1988). However, it would be 
premature to use that correlation for interpreting our observed 
slight trend of helium abundance with luminosity. The metal- 
licity plays an important role both for radiative levitation and 
the radiation-driven wind. It would therefore be of enormous 
interest to determine metalUcities for our high-resolution sdB 
spectra and use them as input for theoretical calculations of ra- 
diative forces and wind mass loss, which, however, is beyond 
the scope of this paper. 
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4.4. Peculiar Ha profiles 

All four stars with peculiar Hor profiles have a higher luminos- 
ity than the majority of sdB stars (see inset in Fig. 9). This gives 
us a hint to a possible explanation: Heber et al. (2003) observed 
similar Ha profiles in the four most luminous objects of their 
sample, and suggested stellar winds as the cause. Similar to our 
data, their profile shapes show broadening around the core re- 
gion, while the line depth is lower than in the model spectrum. 
The core itself is flattened in one case, and shows a faint emis- 
sion in the other cases, consistent with our own observations. 

First calculations of Ha Unes from spherical NLTE mod- 
els including mass loss have recently become available (Vink 
2004). While they can explain the occurrence of small emis- 
sions in the core and, therefore, a weakening of the core, they 
are unable to match the observed widths of the line wings. In 
the latter respect, there is no improvement yet compared to 
the static, plane paraUel NLTE models we use. Therefore the 
physical efl'ects causing the extraordinary line broadening re- 
main obscure. Obviously, there is a need for more sophisticated 
model atmospheres. 

In the sdB mass loss calculations mentioned above, it is 
possible to reproduce the emission with mass loss rates of 
10"' ' M0/yr. On the one hand, such a strong mass loss agrees 
with our conclusion from the previous subsection that mass loss 
rates above 10"'^Mo/yr would be needed to explain the up- 
per part of the heUum abundance range of sdBs. On the other 
hand, the helium abundance of the four stars discussed here 
Ues between -3.41 < \og(y) < -1.60, leading us to recaU the 
statement that steUar winds need not automatically cause a cor- 
relation of luminosity and heUum abundance. 




log(L/L,,,) 

Fig. 10. Cumulative luminosity function: cumulative normalized 
counts versus luminosity in units of the Eddington luminosity. 
Shown are our data (filled circles) along with the measurements of 
Edelmann et al. (2003, open circles), Maxted et al. (2001, open 
triangles), and Saffer et al. (1994, open squares). ZAEHB and 
TAEHB are from Dorman et al. (1993) for metallicities [Fe/H] = 
0.00, -0.47,-1.48, -2.26 from left to right. 



4.5. Comparison with other samples 

SdB samples of comparable size were presented in the last 
decade by Saffer et al. (1994), Maxted et al. (2001) and 
Edelmann et al. (2003). Their main selection effect is the 
magnitude-limitation of the surveys that provided the targets 
(Falomar-Green Survey, HQS). In contrast, the selection effects 
of our sample are not weU known. Since the SPY targets were 
drawn from various catalogs with criteria for selecting potential 
white dwarfs, the hot subdwarfs were included only by mis- 
classification, rendering a definition of sdB selection criteria 
impossible. We now attempt to investigate whether the selec- 
tion effects of our sample are similar to those of the samples 
just mentioned, or whether any systematic differences can be 
found. 

An adequate means for comparing different samples of sdB 
stars is their respective cumulative luminosity function, shown 
in Fig. 10 for our sample along with the three studies named 
above. It can be seen that the SPY data agree quite well with 
the observations from Maxted et al. (2001) and Edelmann et al. 
(2003), whereas there is a larger average offset to the data 
of Saffer et al. (1994), which has already been discussed in 
Edelmann et al. (2003). Furthermore, several of our stars lie be- 
low the ZAEHB, which is plotted in Fig. 10 for various metal- 
licities. Interestingly, the percentage of stars below the solar 
metaUicity ZAEHB (12%) is in better agreement with the ob- 



servations from Saffer et al. (1994, 15%) - despite the overall 
offset between the two functions - than with the other two stud- 
ies (3% and 5%). Similarly, at the highest luminosities, our data 
again agree best with the sample of Saffer et al. (1994). 

The observed moderate differences between the samples 
are most probably the result of differences in resolution, signal- 
to-noise ratio, homogeneity, and models used for line profile 
fitting. As far as the latter are concerned, one must not only dis- 
tinguish between the actual use of LTE or NLTE models (Saffer 
et al. 1994 only use the former), but also between the param- 
eter regions in which they are appUed: Edelmann et al. (2003) 
and Maxted et al. (2001) use NLTE models for objects hav- 
ing Jeff > 27 000 K, our study for sdBs with T^ff > 32 000 K. 
As another example of model differences, we mention that the 
LTE models of Saffer et al. (1994) did not include metal line- 
blanketing, while our LTE models do. Differences between 
NLTE and LTE model atmosphere results may be as large as 
Alog(L) = 0.1 dex (see Heber & Edehnann 2004). 

We conclude from the overall comparison of the samples 
and the lack of excessive deviations that the selection effects of 
our study are mainly similar to those of previous studies - how- 
ever, as far as flux limitation is concerned, we attempt to make 
a more reliable statement below. A more detailed discussion in 
the context of stellar evolution is deferred to Sect. 6. 
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Fig. 11. Cumulative luminosity function for all our objects (filled cir- 
cles, as in Fig. 10), as well as for three magnitude- selected subsets of 
equal number of stars: B < 14.7 (open circles), 14.7 < fi < 15.6 (open 
triangles), B > 15.6 (open squares). 

4.6. Flux limitation 

Our programme stars are drawn from flux-limited samples. 
This may introduce an observational bias against objects of 
lower luminosity, although the limiting magnitudes of the sur- 
veys (e.g. HES and HQS) are considerably deeper than the limit 
for SPY. If such a bias were present it would be more difficult to 
reconcile the lower end of the cumulative luminosity function 
with evolutionary predictions, since we would have to correct 
for low luminosity stars missing in our sample. 

We therefore carried out a test by comparing the cumula- 
tive luminosity functions for three subsets of equal number of 
stars to that of the fuU sample (see Fig. 11). A bright subset 
(B < 14.7), a medium bright subset (14.7 < B < 15.6), and 
a faint subset (B > 15.6) were defined. As can be seen from 
Fig. 11, the cumulative luminosity function of the faint sub- 
set agrees very well with that of the bright and medium bright 
subsample, and all three agree well with that of the entire sam- 
ple. We therefore conclude that no bias due to flux Umitation 
is present. This makes the cumulative luminosity function a 
strong tool for the comparison of observation with predictions 
by the theory of stellar evolution (see Sect. 6). 

5. Double-lined objects 

5.1. Companion classification 

For the 24 sdB stars showing spectral signatures of a cool 
companion, accurate values of the stellar parameters cannot 
be expected from a simple model spectrum fit. In particular, 
the surface gravity is probably underestimated in many cases. 
Nevertheless, we performed Une profile fits of the Unes ¥16 



and bluewards for 20 objects, attempting to minimize the ef- 
fects of companion contribution, since the latter decreases with 
decreasing wavelength. In order to estimate the helium abun- 
dance, we had to include He 1 4471 A in some cases, when no 
other helium line was present in the blue part of the spectrum. 
Our goal was to examine if the fit results can still be useful 
to determine absolute visual brightnesses My sdB of the sdBs, 
from which we can estimate the companion spectral type using 
V and J of the sdB-nmain sequence system. 

For testing purposes, we "contaminated" three of our sdB 
stars showing no companion signatures with a certain amount 
of fight from a cool companion. Since the SPY stars have not 
been flux caUbrated, it was first necessary to create a typi- 
cal spectral energy distribution for the sample stars. This was 
achieved by using flux cafibrated low resolution spectra of sdB 
stars from Salomon (2003), which were chosen to have T^ff- 
and log (g)-values similar to the sample stars. All our spec- 
tra were convolved with a Gaussian of 2.5 A FWHM and then 
rebiimed to 1.0 A in order to match the spectra used for the 
companions (see below). After dividing the sample stars by the 
corresponding flux calibrated ones, a fit to the continuum of 
the resulting spectrum was performed manuaUy. Finally, each 
sample star spectrum was divided by the respective fit curve, 
resulting in a pseudo-flux calibration for them. 

For typical main sequence stars, we used spectra from the 
STELIB^ database (Le Borgne et al. 2003) for companion types 
F9, G4, KO, and K2. Their resolution is A/I < 3 A, which 
matches the treatment of our spectra described above. From 
the absolute visual brightness of the sdB stars and the cho- 
sen companion spectrum, flux contribution and total magni- 
tudes in B, V, and J were calculated (Bessell & Brett 1988; 
Aller et al. 1982; Cox 2000), and composite spectra were 
produced (see Fig. 12). Now, Balmer Une profile fitting for 
H6 and bluewards was performed, keeping the helium abun- 
dance fixed for simplification. My sdB was then calculated from 
the derived values as described in Sect. 4.2. Those values 
show that the efl'ect of companion contribution can be quite 
large for F-type (Alog(^) = -0.23 to -0.6 dex) and G-type 
(Alog(g) = -0.06 to -0.29 dex) companions, which justifies 
our separate treatment of single-lined and double-lined sdB 
stars. For K-type companions the gravity estimate is almost un- 
affected (A log (g) = -0.02 to -0.11 dex). 

Let us now turn to the composite spectrum sdBs from SPY. 
We do not attempt to deconvolve the spectra, but analyze their 
photometric data (B, J). For obtaining apparent brightnesses 
wytot for the system and my sdB for the sdB alone, we begin 
with adopting (B - V)tot = -0^25 ^ to the system, the typical 
value for a single sdB star (Stark & Wade 2003), and assume 
ms tot = otb sdB . After having performed the calculation of com- 
panion type once, as described below, we then correct (B- y)tot 



^ http://webast.ast.obs-mip.fr/stelib/ 

^ Although the Hamburg/ESO Survey provides B - V colours 
(ChristUeb et al. 2001), they have been calibrated only for main se- 
quence stars with B - V > -CTl to an accuracy of ±(^1. Since an 
even better accuracy would be necessary for our purposes, a calibra- 
tion for sdB stars would have to be carried out beforehand, which is 
beyond the scope of this paper. 
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Fig. 12. Observed and artificial composite spectrum. Both spectra 
were scaled to the same flux level, and then displayed with an off- 
set. Obsetyed spectrum: HE 1309-1102, derived companion type G9, 
corrected MvsdB = 4?0. Artificial spectrum: HE 2349-3135 with G4 
companion, MvsdB = 3?86. 

according to the calculated type and start over again with a new 
mvtot- 

The iteration itself is done in the following way: The dis- 
tance modulus of the system is calculated from Mysds and 
"JvsdB- The J magnitude is determined by subtracting the typi- 
cal value of V - J = -0^6 for single sdB stars (Stark & Wade 
2003, consistent with our own findings). Then, by compari- 
son with the measured mjtot, wjcomp is calculated, leading to 
A^jcomp when combined with the distance modulus. This au- 
tomatically yields a companion type when compared to litera- 
ture (Bessell & Brett 1988; Aller et al. 1982; Cox 2000), pro- 
viding us with a value for Mvcomp and subsequently mvcomp- 
Now, a more reliable my sdB was calculated from the latter and 
my lot, and the next iteration step was performed. The iteration 
is stopped when the spectral subtype of the companion stays 
the same from one step to the next, reflecting our maximum 
achievable accuracy. 

The iteration leads to good estimates of the companion 
types for all our artificial composite sdB stars, even when the 
companion contribution to the total flux is quite high. The dif- 
ference of input and output type lies between and 2 subtypes. 
This clearly shows that the method can indeed be appHed to 
our observed sample of composite sdB stars. However, we must 
point out three sources of error: 

1. We adopt a typical value of (V-J)sdB = -07'6. Since the ob- 
served scatter is quite large (see Fig. 2 and Stark & Wade 
2003), this could produce errors in our companion types. 
The estimated error for B, from which V is determined (see 
below), is 0^2, J errors from 2MASS he at about 0^1. For 
A(y - 7)s(iB = ±0'?'2, the difference is two subtypes or less 



for low companion contribution (i.e. K types), and even 
smaller for brighter companion stars (i.e. earlier types). 

2. Similarly, we adopt a typical value of (B - V)sdB = -0^25, 
which in reaUty also shows some small scatter of the order 
of 0"1 (Altmann et al. 2004). 

3. For deriving companion types and absolute visual bright- 
nesses, we use typical values given in the literature, again 
neglecting the natural scatter contained in the true values. 

We therefore conclude that the error limits of our results 
are a few subtypes. We do not account for extinction because 
of the relatively high galactic latitudes of our stars. The nine ob- 
jects which are also in the sample of Altmann et al. (2004, see 
Table 1) have 0.004 < Eg-v ^ 0.035, justifying that approach. 
However, extinction might be important for PG 0258-H184 and 
PG 2122-1-157, where Eb-v > 0.1 (Schlegel et al. 1998). It fol- 
lows that for these objects the sdB is somewhat brighter than 
what is derived by our iteration, and hence the companion is 
slightly fainter, and of later type. 

Our derived companion types for the observed composite 
sdB stars are shown in Table 3. They range from F7, which is 
an early-type hmit in our iteration, to K3 for the stars showing 
Mg I in the spectrum. As expected, aU of the flve companions 
that were only deduced from a contribution at Ha are of rel- 
atively late spectral type (Kl to K7). The derived atmospheric 
parameters for the corresponding sdBs are not affected by com- 
panion spectral contribution, since the latter is negligible for 
the lines H6 and bluewards, where line profile fitting was per- 
formed. All composite sdB stars will be investigated in more 
detail in the future, when we will attempt to deconvolve their 
spectra as described in Aznar Cuadrado & Jeffery (2001, 2002). 

5.2. Main sequence or subgiant companions 

The luminosity class of cool companions to sdBs has been 
discussed extensively. For example, Allard et al. (1994) and 
Jeffery & Pollacco (1998) claimed that the companions to their 
sdB stars are mostly overluminous as compared to main se- 
quence stars. In contrast, Aznar Cuadrado & Jeffery (2001, 
2002) found their companion stars to be mostly main sequence 
objects. Recently, Stark & Wade (2003) again discussed the na- 
ture of companions to sdB stars, pointing out that optical and 
infrared photometry alone is not able to distinguish between 
the main sequence and the subgiant scenario. 

Figure 13 shows the absolute visual brightness distribution 
of our non-composite sdB stars. The values range from 3"6 
to 5?0, with two extreme values at 2™76 and 2™94 from one 
post-EHB and one potential HBB object, respectively. In the 
previous section, our iterative calculations of companion types 
yielded as a byproduct absolute magnitudes of the sdB com- 
ponent (see Table 3), which are also shown in Fig. 13. These 
values are in very good agreement with the ones derived for the 
non-composite sdB stars, and are mostly much fainter than the 
typical values of My « 3 for subgiant stars (Aller et al. 1982; 
also see Reid, http://www-int.stsci.edu/~inr/cmd.html, and ref- 
erences therein). This clearly proves that the majority of our 
companion stars are indeed main sequence stars. 



T. Lisker et al.: Hot subdwarfs from SPY - 1. Atmospheric parameters and cool companions of sdB stars 13 

Table 3. Determination of companion types, the results of which are Usted in the second column. The first 15 objects show Mg i in the spectrum, 
the last five stars only show a continuum contribution at Hor. 



Object 


Type 


WBtot 


mjtot " 


™VsdB 


'^V comp 


(fi - V)to, 




Remark 






mag 


mag 


mag 


mag 


mag 


mag 




HE 1038-2326 


F7 


15.9 


15.22 


16.7 


16.4 


+0.11 


4.0 


F7 or earlier 


HE 1140-0500 


Gl 


14.8 


15.17 


15.1 


17.3 


-0.18 


2.4 




HE 1221-2618 


GO 


14.7 


13.91 


15.4 


15.2 


+0.13 


4.6 




HE 1309-1102 


G9 


16.1 


16.14 


16.4 


18.2 


-0.13 


4.0 




HE 1352-1827 


G3 


16.0 


15.75 


16.5 


17.2 


-0.03 


4.1 




HE 1422-1851 


F9 


16.3 


16.17 


16.7 


17.7 


-0.07 


3.2 




HE 1441-0558 


G5 


14.4 


13.79 


15.0 


15.2 


+0.08 


4.9 




HE 2156-3927 


K3 


14.26 


14.44 


14.5 


17.2 


-0.18 


4.1 


Altmann et al. (2004): (B - V)tot = -0.19 


HE 2322-0617 


G9 


15.7 


15.58 


16.1 


17.4 


-0.09 


4.3 




HE 2322-4559 


G4 


15.5 


15.25 


16.0 


16.8 


-0.03 


4.2 




HS 1536+0944 


KO 


15.6 


15.27 


16.0 


17.2 


-0.07 


4.7 




HS 2216+1833 


Gl 


13.8 


13.32 


14.4 


14.7 


+0.04 


4.2 




PG 0258+184 


G8 


15.3 


14.99 


15.7 


16.7 


-0.06 


4.4 




PG 2122+157 


F7 


15.0 


13.95 


16.4 


15.0 


+0.28 


5.1 


F7 or earlier 


TON S 155 


F7 


16.1 " 


15.69 


16.8 


17.0 


+0.02 


3.5 


F7 or earlier 


HE 1200-0931 


Kl 


16.2 


16.07 


16.5 


18.1 


-0.11 


4.6 




HE 1254-1540 


K7 


15.2 


15.56 


15.5 


19.0 


-0.22 


4.3 




HE 1419-1205 


K5 


16.2 


16.49 


16.5 


19.6 


-0.21 


4.2 




HS 2125+1105 


K4 


16.4 


> 16.5 


16.7 


19.2 


-0.18 


4.5 


K4 or later 


KUV 01542-0710 


K2 


16.3" 


16.30 


16.4 


18.8 


-0.18 


4.0 





" Data from 2MASS, obtained by usage of the VizieR database (Ochsenbein et al. 2000). 
" Johnson V magnitude 
" Stromgren y magnitude 



In addition, our values of My sdB are supported by observa- 
tions of twelve EHB stars in the globular cluster NGC 6752. Its 
known distance allows a direct determination of their absolute 
brightnesses, which lie in the range 3™80 < Mv sdB < 4™66, 
with a median value of My sdB = 4?21 (Moehler et al. 1997). 
The median value for our non-composite objects is My sdB - 
4"20, and My sdB = 4™2 for the sdB component of the compos- 
ites. 

We therefore conclude that the vast majority of cool com- 
panions to observed sdB stars are main sequence stars, and 
consider the question of companion luminosity class to be set- 
tled. However, it follows automatically that sdB+subgiant sys- 
tems can hardly be detected because of the subgiant being 
brighter than the typical sdB. Since "our goal should be to 
describe Nature, not to confirm our own favorite prejudices" 
(R. A. Wade), it has to be pointed out that such systems may 
exist, yet mostly remain undetected. 



5.3. Helium abundance 

Aznar Cuadrado & JefFery (2002) find that the composite- 
spectrum sdB stars show lower photospheric helium abun- 
dances than the non-composite objects, the former having 
log (y) < -2. Our derived helium abundances for the 20 com- 
posite sdB stars for which fitting was performed can be seen 
in Table 1. They are expected to be lower limits, since the 
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Fig. 13. Shaded histogram: Distribution of absolute visual bright- 
ness of non-composite objects. The two leftmost counts come from 
the post-EHB star HE 0415-2417 and the potential HBB star HE 
0151-3919. Open histogram: Distribution of 20 out of 24 sdBs with 
cool companions, as derived from the iteration described in Sect. 5.1. 
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Fig. 14. Shaded histogram: Distribution of photospheric helium abun- 
dance of non-composite sdB stars, as in Fig. 7. Open histogram: 
Distribution of 20 out of 24 objects with cool companions. 

presence of a cool companion contributing to the total flux di- 
lutes the helium lines and therefore decreases the helium abun- 
dance derived from spectral fitting. A similar test as described 
in Sect. 5.1, using HE 2237-1-0 150, shows that the helium abun- 
dance is underestimated by 0.17 dex for an F9 companion and 
0.03 dex for a K2 companion. 

Figure 14 contains the derived values as a histogram to- 
gether with the results for non-composite sdB stars. It is clear 
that the helium abundances for composite-spectrum sdB stars 
cover the whole range of abundances from log {y) < -4 to 
log iy) > -1, witii 50% of the stars having log (y) > -2, in ob- 
vious disagreement with the observations of Aznar Cuadrado 
& Jeff^ery (2002). 

The stars discussed here have formed in binaries with rel- 
atively large orbital separations and long periods (1st RLOF 
channel, see Sect. 6). Aznar Cuadrado & Jeffery (2002) sug- 
gested that the different strengths of tidal effects in short-period 
and long-period sdB binaries may cause differences in the at- 
mospheric helium content of these two groups. This argument 
was taken up by Edelmann et al. (2003) as a possible expla- 
nation for their two distinct sdB sequences in the Teff-log (>)- 
plane (see Fig. 8). The absence of any differences between 
single-lined and double-lined sdBs in our sample rules out 
such a correlation. Furthermore, a large number of the non- 
composite sdB stars may be single objects, which is not taken 
into consideration in these discussions. 

6. Observation versus theory 

In order to draw conclusions about potential formation scenar- 
ios and the evolutionary status of our objects, we consider two 
theoretical studies. Dorman et al. (1993) computed evolution- 



ary tracks for the canonical core mass established by the helium 
flash in the core of a red giant star. The core mass depends on 
metallicity, since the elemental abundances determine the core 
mass necessary for helium ignition on the first giant branch 
(FGB). Dorman et al. (1993) used various metallicities and en- 
velope masses, defining the standard ZAEHB and TAEHB, as 
well as EHB and post-EHB evolution, which we have already 
used in Sect. 4.5. We will take those calculations as representa- 
tive for scenarios of single star sdB evolution, which we define 
in such a way that any companion star that might be present 
does not affect the evolution of the sdB. 

In contrast, the binary population synthesis calculations of 
HPMM aim to distinguish between several binary formation 
channels of sdB stars and then simulate the evolution of the 
resulting subsamples. Those calculations are the most recent 
and most extensive simulations for binary formation scenarios 
of the sdB population. They can be divided into three channels: 

(i) The common envelope ejection channel. Dynamically un- 
stable mass transfer in a binary system with one star being 
on the FGB results in the formation of a CE. Subsequently, 
the spiraling-in of both stars and finally the ejection of the 
CE can lead to an sdB in a very close binary system, if the 
giant's core is still able to ignite helium. The second mem- 
ber of the system can either be a main sequence star ("1st 
CE ejection channel") or already a white dwarf ("2nd CE 
ejection channel"). 

(ii) The stable Roche lobe overflow channel. This channel is 
analogous to the CE ejection channel, except for the mass 
transfer being stable. The value of ^crit decides which of the 
two channels applies to a system: only if the ratio of red gi- 
ant mass to companion mass is below q„[i (adopted to be 
either 1.2 or 1.5 in the simulations) is mass transfer stable. 
This results in binary systems with a much larger separa- 
tion and therefore much longer periods, since there is no 
spiraling-in phase affecting the orbital parameters. Again, 
one has to distinguish between the "1st RLOF channel" and 
the "2nd RLOF channel". However, as shown in HPMM, 
the 2nd RLOF channel (stable mass transfer onto a white 
dwarf) is negligible because white dwarfs of sufficiently 
high masses are extremely rare. 

(iii) The merger channel. The merger of two He-WDs caused 
by loss of energy due to gravitational wave radiation can 
produce an sdB star, if the mass is sufficient for core helium 
burning. Unlike sdB stars from other chaimels, the objects 
produced here are single stars. 

To examine the effects of various poorly known physical 
parameters in the scenarios just described, HPMM produced 
twelve simulation sets for different values of the following 
quantities: CE ejection efficiency acE, fraction ath of CE ther- 
mal energy used for its ejection, critical mass ratio ^crft above 
which mass transfer is stable on the RGB, initial mass ratio 
distribution of the progenitor binary systems, and metallicity. 

HPMM compared the periods and minimum companion 
masses for the single-lined sdB binaries observed by Maxted 
et al. (2001) and Morales-Rueda et al. (2003) to the predictions 
of their numerical models and selected a simulation set that 
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Fig. 15. Teff and log (g) values of our sdBs as shown in Fig. 6, along 
with evolutionary tracks calculated by Dorman et al. (1993). Tracks 
are shown for solar metallicity (heUum core mass Mcore = 0.469 M©, 
light grey lines) as well as for [Fe/H] = -2.26 (Mcore = 0.495 M©, 
dark grey lines). The respective envelope masses M^nv are displayed in 
the figure. The solid part of the tracks marks the core helium burning 
phase, while this has already ceased when the star is on the dotted part 
of the tracks. 
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Fig. 16. Cumulative luminosity function of our sdB stars as shown 
in Fig. 10, along with sdB luminosity evolution with time as calcu- 
lated by Dorman et al. (1993). The evolutionary tracks (grey lines) are 
followed from the ZAEHB until they leave our Teff-log (^)-parameter 
space (T^ff < 40 000 K, log (g) > 4.5). For these tracks, the metallici- 
ties are from left to right: [Fe/H] = 0.00, -0.47, -1.48, -2.26 (corre- 
sponding to core masses of Mcore = 0.469, 0.475, 0.485, 0.495 M©). 
For the leftmost track the envelope mass is Mgny = 0.002 M©, for the 
other tracks Menv = 0.003 M©. 



matches those observations best ("best-fit model"). However, 
this comparison does not include stars which formed via the 
merger channel, since these are single objects. Hence, it is very 
interesting to compare HPMM's predictions directly with our 
spectroscopic results, because our sample is not biased against 
single stars. 

To make a comparison possible and reasonable, it is im- 
perative to consider observational selection effects and apply 
them in the same way to the theoretical sample. HPMM showed 
that the so-called GK selection effect is the most important 
one, since observational studies typically select against G, K, 
and earher type main sequence companions. Those either out- 
shine the sdB or cause a composite spectrum which renders the 
analysis of the sdB difficult or even impossible. In the simula- 
tions, HPMM apply this effect by excluding all sdB +main se- 
quence systems with the companion either being brighter than 
the sdB or having T^g > 4000 K, corresponding to M and very 
late K types. Since we excluded all composite spectrum stars 
carefully, our sample matches this approach, and allows a de- 
tailed comparison of our data with the theoretical predictions 
of HPMM. 

6.1. Single Star evolution 

Figure 15 compares our results to several evolutionary tracks 
calculated by Dorman et al. (1993) for different core and en- 
velope masses. It can be clearly seen that at least one of our 



stars (HE 0415-2417 at T^g = 32768 K, log(g) = 5.12) must 
be in the post-EHB stage of evolution (dotted part of the tracks 
in Fig. 15). For a more thorough comparison with these the- 
oretical predictions, we plot in Fig. 16 our cumulative lumi- 
nosity function from Sect. 4.5 along with the luminosity evo- 
lution with time given by the described tracks, until they reach 
an effective temperature of 40 000 K. This limit follows from 
our selection criteria. The function should reflect the "speed of 
evolution" across the EHB and in the post-EHB stage. 

It has been argued previously that most stars would be ex- 
pected to he near the ZAEHB (e.g. Edelmann et al. 2003). 
The theoretical tracks, however, show a linear time-luminosity- 
relation while the star is in the EHB strip. The respective 
ZAEHB is defined by the starting point of each track, the 
TAEHB is clearly marked by the sharp bend in the tracks. We 
therefore have to expect a homogeneous distribution of sdB 
stars through the EHB strip, if there are no observational se- 
lection effects. 

While the discussed evolutionary tracks lie at somewhat 
lower average luminosity, their slope matches our data well, 
except for the six stars at luminosities lower than any of the 
plotted tracks. As shown in Sect. 4.6, any potential effects of 
magnitude limitation are negligible, and do not affect this com- 
parison. Furthermore, replacing NLTE models by metal line- 
blanketed LTE models for the hottest stars would primarily 
lower the luminosity mainly at the top part of the cumulative 



16 



T. Lisker et al.: Hot subdwarfs from SPY - 1. Atmospheric parameters and cool companions of sdB stars 



luminosity function by about 0.06 dex. Adopting a very low 

metallicity would increase the luminosities by 0.03 dex only. 
Hence these systematic errors are unimportant (see Sect. 4.1). 

6.2. Close binary evolution 

We proceed with a detailed comparison of our measurements 
with all twelve simulation sets of HPMM. Our goal is to ex- 
clude several sets and their respective parameter configura- 
tions, to test HPMM's choice of the best-fit model, and to check 
whether these calculations can well reproduce the observed 
sdB distribution. 

A short overview of the input parameters that characterize 
the simulation sets is given now. Three sets in a row (1-1-2-1-3, 
4H-5-I-6, etc.) constitute a sequence in both the common enve- 
lope ejection efficiency acE and the thermal energy fraction ath 
used for ejection. The first set always has qtce = ^th = 0.5, the 
second acE = <^th - 0.75, the third acE - Q'th - 1-0. Sets 1 to 6 
have a critical mass ratio qcnt = 1.5, whereas qan = 1.2 is used 
for sets 7 to 12, allowing stable mass transfer only for a nar- 
rower mass range than the former value. For sets 4, 5, and 6, an 
uncorrelated mass distribution had been adopted for the stars 
in sdB progenitor binaries, whereas in all other sets, HPMM 
used a constant distribution of the progenitor system mass ra- 
tio. Sets 10, 1 1, and 12 were calculated with a typical thick disk 
metallicity (Z = 0.004) while the other sets are for population I 
metallicity (Z = 0.02). 

6.2.1. The reff-log(g)-plane: visual inspection 

Figure 17 shows the simulated data from HPMM in the T^g- 
log (§)-plane. The number density is grey-scale coded; higher 
densities of sdB stars correspond to darker grey shading. Note 
that the total number of simulated stars is not the same for 
different simulation sets. Overlaid on each simulation data set 
are our measurements. The simulations shown here include the 
GK selection effect, as explained above, and are therefore well 
suited for comparison with our observations. 

Sets 4, 5, and 6 do not match our data, since their pre- 
dicted sdB density at higher temperatures is much lower than 
for cooler stars, in obvious contrast to our measurements. Sets 
3, 9, and 12 all show a significant number of objects at low T^g 
and higher log (g), i.e. in the lower right area of the diagram. 
These stars have somewhat lower luminosities, and are not seen 
in our data. Since no bias due to flux limitation is present in our 
sample (see Sect. 4.6), we can state a trend against these simu- 
lation sets. 

From visual judgement only, sets 1, 2, 7, 8, 10, and 11 all 
match the observed sdB distribution well. In all sets, a trend can 
be seen that the average surface gravity at lower temperatures 
is somewhat higher in the simulations than in our data. Again, 
this cannot be explained with a potential flux limitation bias, 
therefore, the observed shift could point towards the necessity 
of refining some model details. 

Two stars are not matched by any simulation set: HE 
0151-3919 at T^b = 20841 K, log(^) = 4.83 and HE 
0415-2417 at T^b = 32768 K, log(^) = 5.12. We cannot rule 



Table 4. Likelihood ranking of the HPMM simulation sets. The log- 
arithmic ratio of the likelihood for the best matching set 10 to the 
likelihood for each set is given. 
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out that the former is a blue HB star, which will ascend the 
AGB after core helium exhaustion and therefore may belong to 
a different evolutionary group than the EHB stars. The latter is 
believed to be in the post-EHB stage (see Fig. 15), which lasts 
only a fraction of the time spent on the EHB. Since there are 
no additional stars observed in this parameter region, we do not 
consider this to be a discrepancy with theory. 

In Sect. 4.3 we found that at higher Teff the sdB stars are 
somewhat more luminous, and tend to have a higher helium 
abundance. The corresponding T^s-log (^)-range of simulated 
data points is dominated by stars that formed in the merger 
channel. Since most of the hydrogen is expected to be burned 
up during the merging process, those objects are more likely 
to have a higher helium abundance than sdB stars from other 
channels (Han et al. 2002). This qualitative agreement with 
our results can be tested in a quantitative way by using ra- 
dial velocity variations to distinguish between single objects 
from the merger channel and binaries (Napiwotzki et al. 2004, 
Napiwotzki et al. in prep.). 

6.2.2. Tlie TeflF-log (g) plane: statistical analysis 

To check the visual judgement, we performed a statistical test 
for the match of the different sets, using the maximum like- 
lihood method. Thereby, we adopt as probability for each of 
our sdBs the fraction of simulated stars in a box around the 
data point, and then multiply the probabilities of all points. The 
size of the box corresponds to two times our errors for stars 
with one exposure, which was carefully chosen to avoid ran- 
dom fluctuations on the one hand, and too much smoothing on 
the other hand. In Fig. 17 the bin size is only half of that value 
- i.e. it equals the errors - and one can clearly see that quite 
large fluctuations can occur, since the total number of simu- 
lated stars only lies between 1000 and 2000, and is distributed 
over several hundred bins. 

The ranking derived from the likelihood of each set to 
produce the observed distribution is given in Table 4, and 
clearly favors set 10, which is characterized by a low efficiency 
(orcE = ath = 0.5), low metallicity (Z = 0.004), and a con- 
stant mass ratio of the progenitor binaries. We want to point 
out that the given likelihood ratios cannot be exact, since the 
simulation data are provided as discrete number counts rather 
than a smooth distribution function. Nevertheless, the overall 
ranking can be considered reliable. Sets 8, 11, and 2 form a 
group that still matches the observations well, while all of the 
following sets give a poor match to our data. Note that HPMM 
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Fig. 17. Comparison of our derived values for T^e and log (g) (symbols and typical errors as in Fig. 6) with the simulations from HPMM. The 
theoretical predictions are shown as shaded Teff-log (g)-boxes, where a higher sdB density per box corresponds to darker shading. The grey 
scale is shown below the figures. Continued on next page. 



favored sets 2 and 8, where the former is their best-fit model. The rejection of sets 4, 5, and 6 is nicely confirmed, since they 
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rank last, together with set 9. Similarly, the trend that we stated 6.2.3. The Cumulative luminosity function 

above against sets 3, 9, and 12 is supported by the statistical 

analysis. p.^ compares the cumulative luminosity function from all 

simulation sets of HPMM with our data. It is obvious that in all 
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Fig. 18. Cumulative luminosity function of our sdB stars as shown in 
Fig. 10, along with the functions given by the HPMM simulation sets. 

panels the respective leftmost track gives the worst match to the 
data. These tracks correspond to sets 3, 6, 9, and 12, which we 
have already found a trend against. Their common parameter 
values are qtce = »th - 1 , leading to the conclusion that CE 
ejection processes do not reach 100% efficiency. 

In contrast to the luminosity functions derived from the 
Dorman et al. (1993) models, the best matching simulation sets 
of HPMM - as derived from our statistical analysis - are con- 
sistent with the observations at the low luminosity end of the 
luminosity function. It can be clearly seen, however, that they 
predict a smaller slope than we observe, and that even set 10 
is far from reproducing the observed function. As mentioned 
before, any systematic errors due to LTE/NLTE or metalUcity 
differences would not affect this statement. 

6.3. Other sdB samples 

We also investigated the sdB samples studied by Saffer et al. 
(1994) and Edelmann et al. (2003). A comparison to the theo- 
retical luminosity function in the context of the single star sce- 
nario (as described for our sample in Sect. 6.1) can be made 
from Figs. 10 and 16. It qualitatively shows the same overall 
behaviour as we found for our survey, i.e. the slope of the cu- 
mulative luminosity function is the same, but the offsets of the 
luminosity scale are different, and the numbers of stars of too 



low or too high luminosity differ. We also compared these sdB 
samples to the HPMM simulation sets by visual inspection of 
the TefF-log (^)-plane (as described in Sect. 6.2.1) and the cu- 
mulative luminosity function (as described in Sect. 6.2.3). In 
general the result is the same as for the SPY sample. The lu- 
minosity function favors models with uncorrelated mass ratios 
(in particular set 4), while the Teff-log (^)-plane rules them out. 
The latter favours the same models (10, 8, 2, 11) as for SPY. 
This directly leads us to a discussion of the current limits of 
sdB formation theories, and to suggestions for future improve- 
ments. 

7. Discussion 

To constrain the two CE ejection parameters further does not 
seem possible nor reasonable because of poor knowledge of 
CE ejection physics. This is also reflected in the controversial 
statements by HPMM and Soker & Harpaz (2003) about the 
treatment of CE ejection processes. The latter authors criticize 
HPMM's criterion for CE ejection, and suggest to account for 
enhanced mass-loss rates by binary interaction, rather than in- 
cluding CE ionization energy as HPMM did. 

To decide between a value of 1 .2 and 1 .5 for q^nt would 
be of enormous interest for the predicted number of sdB stars 
that are hidden from view because they are outshone by main 
sequence companions of spectral types B to F. Especially the 
percentage of A-type main sequence stars in the Galaxy that 
have sdB companions is highly sensitive to ^crit> predicted to be 
0.75%if^crit = 1.5, but only 0.19% if 9„it = 1.2. HPMM point 
out that "the effects of a tidally enhanced wind can to some de- 
gree be impUcitly included by using larger values of ^crit"- This 
immediately suggests that various values of q^h may be real- 
ized in sdB formation, depending on the presence and strength 
of tidal wind enhancement. 

Similarly, it cannot be doubted that various metallicities 
may play a role among the observed sdB stars, since their pro- 
genitors' ZAMS masses range from Mzams < 1 - corre- 
sponding to main sequence lifetimes comparable to the age of 
the Galaxy - to almost Mzams = 3 M0, having a much shorter 
main sequence Ufetime. In addition, sdB stars exist in some 
halo globular clusters, and kinematic studies show that sdB 
stars are present in the thin disk, the thick disk, and the halo 
(Altmann et al. 2004). Hence, a mixture of different parameter 
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configurations and population memberships may be realized in 
nature instead of having fixed values for all physical quantities 
involved. 

These considerations require sdB formation theories to be 
even more complex than the HPMM simulations. In addition, 
the comparison of our data with the predicted cumulative lumi- 
nosity functions of HPMM and Dorman et al. (1993) in the pre- 
vious section can well be interpreted with the necessity of com- 
bining single star and binary formation channels to reproduce 
the observed sdB population. Although not compelling, one 
can see in Fig. 18 that the best-matching slope of the HPMM 
simulation sets is given by set 4 - regardless of the offset - 
which we had immediately excluded because it does not match 
the observed Teff-log (^)-distribution. Since the slope of the cu- 
mulative luminosity function in the single star scenario (see 
Fig. 16) is steeper, we could remedy the discrepancy in the bi- 
nary scenario (HPMM) by allowing for a contribution by a sin- 
gle star evolution channel, i.e. both scenarios contribute to the 
observed sample. An alternative possibility is that our sample 
is biased against stars of higher luminosity, which, if added to 
the function, would flatten the steep slope and probably smooth 
away the sharp bend at the high luminosity end. Although this 
seems unlikely for our sdB sample, we still have to investigate 
the relevance of the sdO stars in the SPY survey for the EHB 
evolution, since a close connection of the hydrogen-rich sdO 
stars to the EHB has often been claimed. Note that sdO stars 
have also been excluded in the investigations of other pubhshed 
samples. 



8. Conclusions 

We have presented the results of a spectral analysis of 76 sdB 
stars from SPY. 24 objects show spectral signatures of a cool 
companion, which we investigated from optical and infrared 
photometry. The majority, if not all, of double-lined stars have 
main sequence companions of types F to K. When focusing 
on helium abundances of composite sdB stars, we find no dif- 
ference when compared to non-composite objects. Of the 52 
single-lined stars, four show pecuUar Ha profiles, possibly in- 
dicating stellar winds. The luminosity distribution of our sam- 
ple is found to be in good agreement with previous studies. 
The tracks for luminosity evolution of single stars calculated 
by Dorman et al. (1993) agree in shape with the cumulative 
luminosity function of our stars on and above the EHB. They 
show, however, a sUght off'set in luminosity, and in addition 
they cannot explain the objects below the ZAEHB. The best- 
matching simulation sets of the binary population synthesis cal- 
culations by HPMM match the observed sdB distribution in T^g 
and log (g) very well, but they cannot reproduce the cumulative 
luminosity function of our stars. The mismatch is most obvi- 
ous for the simulation sets using 100% efficiency of CE ejec- 
tion (sets 3, 9, and 12), which also holds for a comparison with 
the luminosity function of other observational datasets (Saff'er 
et al. 1994; Edelmann et al. 2003). Furthermore, there is some 
evidence against an uncorrelated mass distribution of the pro- 
genitor systems, which follows from the trend of rejecting sets 
4, 5, and 6. 



We conclude that a combination of single star and binary 
formation channels would be necessary to achieve full under- 
standing of sdB formation processes. In order to solve the prob- 
lems that the latter still pose, future simulations are required 
to incorporate a more sophisticated description of the physics 
involved. In parallel, high-quality observational samples hke 
SPY are necessary to enable an even better judgement of the 
simulated sdB population. A crucial test of theoretical calcula- 
tions would be the comparison of predicted and observed frac- 
tion of radial velocity variables, i.e. stars that formed in the CE 
ejection channels and thus are close binaries. Since SPY was 
initiated to search for RV- variations, the data are well suited for 
such an investigation. The subsequent step will then be to de- 
termine orbital parameters of the RV-variable sdBs from SPY 
by measuring their radial velocity curves. First results of these 
projects are reported by Napiwotzki et al. (2004), and more de- 
tails will be presented in Napiwotzki et al. (in prep.). 

Equally important is the completion of the analysis of hot 
subdwarfs from SPY, i.e. the investigation of the various sdO 
stars observed. This class is much less homogeneous than the 
sdBs, and it remains to be seen whether they (or at least some 
of them) are intimately connected with the EHB. Increasing the 
sample size of sdB and sdO stars should then be the ultimate 
step. For a detailed comparison with evolutionary models, it is 
mandatory to improve number statistics. A huge number of sdB 
and sdO stars is now being discovered by sky surveys such as 
the Sloan Digital Sky Survey, and needs to be studied. 
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